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It has recently been shown that light can be stored in Bose-Einstein condensates for over a second. 
Here we propose a method for realizing a controlled phase gate between two stored photons. The 
photons are both stored in the ground state of the effective trapping potential inside the condensate. 
The collision-induced interaction is enhanced by adiabatically increasing the trapping frequency and 
by using a Feshbach resonance. A controlled phase shift of n can be achieved in one second. 
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Photons are ideal carriers of quantum information over 
long distances. It is interesting to explore their poten- 
tial for the implementation of quantum information pro- 
cessing as well. This is particularly relevant for quan- 
tum repeaters [1-3], which would allow one to distribute 
quantum states over distances that are inaccessible by 
direct transmission. Quantum repeaters require both 
the capability to store photons for relatively long times 
and to perform efficient quantum gates between them 
0. Potential architectures where storage and quantum 
gates can be achieved in the same system are particu- 
larly attractive. Recently it was shown that light can be 
stored for over a second in a Bose-Einstein condensate 
(BEC) 0, making condensates a very interesting can- 
didate system for the implementation of quantum mem- 
ories. Quantum repeaters can tolerate long gate times 
in the sub-second range, since repetition rates are in 
any case limited by other factors such as communication 
times and transmission probabilities. It is therefore of 
great interest to explore the potential for photon-photon 
gates in BECs, where interactions between stored exci- 
tations are weak, but non-zero. 

In the following we describe a concrete proposal for 
realizing such photon-photon gates in BECs. Our work 
builds on Refs. 0,0], but we focus on the case of two sin- 
gle photons interacting. In this extremely low-intensity 
regime achieving significant controlled phase shifts is not 
straightforward. However, we show that phase shifts of n 
can be achieved on sub-second timescales by combining a 
Feshbach enhancement of the relevant scattering length 
and an adiabatic compression of the trap after the light 
has been stored. The fidelity of photon-photon gates can 
be affected by unwanted multi-mode effects, see e.g. Rcf. 
[3] ■ In the present proposal these effects are greatly sup- 
pressed by the fact that the interaction is much weaker 
than the confinement, ensuring high-fidelity operations. 

Let us assume that the two photons have orthogonal 
polarization. Their propagation inside the BEC can be 
controlled by two independent control beams, leading to 
storage in two different atomic levels 1 and 2, where the 
BEC was prepared in level 0, see Fig. [1] Slow and 
stoppe d lig ht in BECs has been thoroughly investigated 
@, Iil4l0|. Due to the linearity of the equations of mo- 
tion, the physics of storage and retrieval is the same at 




FIG. 1: Level scheme for photon-photon gate. The BEC is 
prepared in level 0. The single photons in modes Ei and 
E2 can be independently stored as delocalized excitations in 
levels 1 and 2, using the control beams Qi and fig. 



the sing le-photon level as for weak classical probe pulses 
111 Il2j . Inside the medium and in the presence of the 
appropriate control beam, the photon is converted into a 
slowly moving polariton, which can be stopped by adia- 
batically switching off the control beam, thus converting 
the photon into a stored atomic spin wave. Running the 
process in reverse allows the reconversion of spin waves 
into photons. Here we focus on the interaction between 
the two spin waves, once the control beams have been 
turned off. Due to the weakness of the collision-induced 
interactions the timescale for the storage and retrieval 
processes is much shorter than the timescale on which sig- 
nificant interaction occurs in the photon-photon regime. 

The dynamics of the atomic spin waves is governed 
by the collisional interactions between atoms in combi- 
nation with the external trapping potential. Spin waves 
in levels 1 and 2 experience an effective trapping poten- 
tial and effective collisional interactions that depend on 
the differences between the atomic scattering lengths in 
the various atomic levels jy]. These differences, which 
are usually small, can be enhanced by Feshbach reso- 
13Ml5|. We consider a situation where both spin 



waves experience the same effective trapping potential, 
and where they are both in its ground state. The lat- 
ter condition can be achieved by carefully matching the 
pulse duration and width of the incoming photons and 
the intensity of the control beams (which determines the 
propagation speed and thus the longitudinal extent of the 
polaritons inside the condensate) to the parameters of 
the effective trapping potential. We focus on the regime 
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FIG. 2: Principle of photon-photon gate, (a) Photons 1 
and 2 are independently absorbed by the BEC. (b) Both 
are converted into atomic excitations that are in the ground 
state of the effective trapping potential, see text, (c) The 
collision-induced interaction between the atomic spin waves 
is enhanced by adiabatically increasing the trapping poten- 
tial, thus reducing the size of the ground state wave functions 
(and of the BEC). (d) The trapping potential is adiabatically 
brought back to its original value, (e) The photons can be 
read out independently. 



where the stored spin waves are localized well inside the 
condensate, cf. Fig. 2. 

The interaction strength, and thus the accumulated 
controlled phase shift for a given time, then strongly de- 
pends not only on the scattering lengths, but also on 
the size of the ground state wave packets. During stor- 
age and retrieval, this size has to be significantly larger 
than a wavelength, due to focusing restrictions for the 
transverse dimensions and in order to justify the slowly 
varying envelope description which underlies the polari- 
ton picture for the longitudinal dimension. However, in 
between storage and retrieval it is possible to adiabati- 
cally increase the trapping potential, thus reducing the 
size of the ground state wave packets while keeping the 
spin waves in the ground state, see Fig. [2l This enhances 
the interaction strength, making controlled phase shifts 
of 7r achievable on one-second timescales. Note that the 
basic ingredients of the present proposal are similar to 
those of single-atom quantum gate s schemes based on 
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cold collisions such as Refs. 



We now describe our proposal in more quantitative 
terms. Our treatment of the spin waves inside the BEC 

The Gross-Pitaevskii (GP) 



18] 



is based on Refs. 
equation for the macroscopic wave- function tpo of the con- 
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where m is the atomic mass, V is the trapping poten- 
tial, Uoo,Uqi,Uo2 are the collisional interaction poten- 
tials, which are related to the corresponding scattering 

lengths 000,001,002 by U oj = n m ° oi , and i/>i,ip2 arc the 
macroscopic wave functions for levels 1 and 2. We will 
make the transition to a single-quantum description for 
the latter in a moment. 

For a sufficiently large condensate, and keeping in mind 
that the perturbation due to the spin waves in levels 1 
and 2 is extremely weak in our case, the solution for tpo 
will be essentially stationary, and the stationary equation 
for a chemical potential /i can be solved in the Thomas- 
Fermi approximation (i.e. neglecting the kinetic term) 
EH, giving 
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{n-V-Uo^f + UwlH 2 ), (2) 



where /i is the chemical potential. This solution of Eq. 
@ can now be inserted into the GP equations for f/h 
and ip2- Corrections to the Thomas- Fermi approximation 
mainly affect the boundary layer of the condensate E3 . 
We therefore expect the present treatment to be correct 
under the above-mentioned condition that the spin waves 
in levels 1 and 2 are localized well inside the BEC. 

In order to describe the few-excitation regime, we re- 
place the macroscopic wave functions tpi , if>2 by quantum 
field operators ip\ , ^2 satisfying commutation relations 
[^i(x), ^t(x')] = 6ij6( 3 \x — x'), in analogy to the tran- 
sition from classical to quantum non-linear optics [20j |. 
They fulfill the equations (neglecting a constant energy 
shift that depends on /i) 

at 2m 
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V" 2 (x) + U22W2 + fWiV'i) V>2,(3) 



where V$ = (1 — f^)V are the effective trapping poten- 
tials and Ua — inh (an — a °' a ° J ) are the effective inter- 
action potentials, which are all modified due to the inter- 
action with the background condensate. These equations 
are analogous to those obtained in Ref. 6] for the two- 
level case. Here we have assumed that the bare trapping 
potential V is the same for all atomic levels. Moreover 
for simplicity we will assume that aoi = &02 implying 
Vi = V2 = V. We require aoi < aoo in order for V to be 
attractive, provided that V is attractive (2ll |. 

Eq. ((3]) allows one to describe the dynamics of any 
number of quantum excitations in levels 1 and 2. How- 
ever, we are interested in the case where there is exactly 
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one excitation in each level. It is then convenient to 
introduce the two-particle wave-function ^^(xi, X2) = 
(O|V>i(xi)-02(x2)| c E > ), where |0) is the state without any 
excitations (i.e. the state where there is just the conden- 
sate in level 0), and 

|$) = /d 3 x 1 d 3 a; 2( /)o(x 1 )^ (x 2 )^ 1 t (x 1 )^(x 2 )|0) (4) 



is the initial state (after storage), which consists of one 
atomic excitation in each level (1 and 2), both of which 
are in the ground state 4>o of the effective trapping po- 
tential V. In the Heisenberg picture for the quantum 
field theory the state remains constant, but the field op- 
erators evolve according to Eq. (|3]). As a consequence, 
the two-particle wave function ^12 defined above evolves 
according to 

iftj^i2(xi,x 2 ,t) = (-^(v? + v 2 ) + V( Xl ) + y(x 2 ) 

+U 12 S^ (x a - x 2 ))Vi 2 (xi, x 2 , t),(5) 

with the initial condition ^12 (xi, x 2 , 0) = 0o(xi)(^o(x 2 ). 
We assume a spherically symmetric harmonic potential 
V(x) = |mw¥, implying </> (x) = (2^) 2 e- 22 ^ . 

It is convenient to transform to center-of-mass and rel- 
ative coordinates defined by X = Xl ^* 2 and r = Xl ^ 2 . 
In these coordinates the wave function is separable at all 



times, -012 (X, r, t) 



! -Vo(X)^(r,i). The center of 



mass wave function exactly remains in the ground state 
of V. The relative coordinate wave function ?/>(r, t) fulfills 
the equation 

i^(r,t) = (-|^V 2 + V(r) + U 12 S (3) (rmr,t), (6) 
ot 2m 

where U12 = t/i 2 2~l. The interaction between the two 
spin wave excitations inside the BEC is thus reduced to a 
fairly simple problem in one-particle quantum mechanics. 

In practice the interaction energy associated with the 
U12 term is two to three orders of magnitude smaller 
than the harmonic oscillator energy scale Tllu. As a con- 
sequence, the use of perturbation theory is well justi- 
fied. Due to the large separation between the two en- 
ergy scales, ip(r, t) remains essentially proportional to the 
ground state, see below. However, there is an energy shift 
relative to the ground state energy, which is given by 

AE = (</>o|^i2<5 (3) (r))l^o) = C7 12 |0 O (O)| 2 = U 12S -\ (7) 
where s = \ is the characteristic length scale of the 

y mu ° 

ground state wave function, which is related to its full 
width at half maximum I by s = J g^ i. This energy 
shift is the basis of our quantum gate proposal. Since 
it is due to the interaction, it only occurs if there are 
two excitations in the condensate, allowing one to real- 
ize a controlled phase gate. The gate naturally has high 
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FIG. 3: Qualitative time dependence of the trap frequency 
uj and of the acquired phase <f). The total time is 2t a +tf, 
where t a is the time required for adiabatically changing the 
trap frequency from ujo to wj or vice versa, and tf is the time 
during which the trap frequency is held fixed at the high value 
wi. The corresponding phases are cj>a and <f>f, giving a total 
phase 2<j> a + 4>f, where </>/ 2> <^ a in the discussed regime. 



fidelity because the correction terms to the ground 
state wave function have amplitudes of order j£ ^ 212. ; 
which is smaller than 10 -2 even for the largest scattering 
length and smallest ground state size that we will con- 
sider. This means that, apart from the phase, the over- 
lap with the initial state remains extremely high, which 
is exactly what is required for high-fidelity operation [3] ■ 
The remaining challenge is therefore to achieve a con- 
trolled phase shift of ir. Let us begin by choosing param- 
eter values that should be straightforwardly achievable. 
For example, one can choose level in the F = 1 sub- 
manifold of 87 Rb, and levels 1 and 2 in the F — 2 sub- 
manifold, giving aoo = 5.39 nm, aoi = 0,02 — 5.24 nm and 
a\2 = 5.58 nm [18j . and a full width at half maximum for 
the ground state wave packet / = 8/xm (corresponding to 
about ten wavelengths) . With these values one finds that 
the time required for a phase of n is 6 minutes, which at 
first sight may seem rather discouraging. We now discuss 
how to overcome this difficulty by acting both on the U12 
factor and the s -3 factor in Eq. ([7]). 



The factor J7 12 = 



V2irh 2 



(a 



12 



) is very small 



for the values given above because there is a quasi- 
cancelation between the two terms in parentheses be- 
cause all the scattering lengths are so similar. A moder- 
ate increase in 012, which can be achieved using a Fesh- 
bach resonance can therefore lead to a very large 

increase of J7i 2 . For example, increasing a\ 2 by a factor 
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of F — 3, which was already demonstrated in Ref. 
for 87 Rb, increases XJ\% by a factor of 24. 

A comparable gain can be achieved by acting on the 
second factor in Eq. ([7]), i.e. on the size of the wave 
function, or equivalently on the trapping frequency. We 
already mentioned in the introduction that I (and thus s) 
cannot be too small during the storage process, because 
focusing becomes too difficult and the slowly varying en- 
velope approximation breaks down. However, the trap- 
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ping frequency can be increased once the photons have 
been stored, see Fig. 2, with the caveat that this increase 
has to be done adiabatically so that the spin waves re- 
main in the ground state of the effective trapping poten- 
tial. The mentioned I = 8/im corresponds to an effective 
frequency Q — 2n 10 Hz, which corresponds to a real trap 
frequency lj — 2?t 50 Hz. This gives a condensate size of 
17 /im for TV = 10 atoms in the Thomas-Fermi approx- 
imation 22j. The effective frequency can be increased 
to Cj = 27r80 Hz in t a = 0.14 seconds while exciting the 
system out of the ground state with a probability that is 
smaller than 0.002. At this frequency the ground state 
size I is 2.9^m and the size of the condensate is 7.4 /mxi. 
For a Feshbach factor F — 3 a phase of order tt can 
then be achieved with tf — 0.73 seconds. Taking into 
account that one has to decrease the frequency before 
readout, the total gate time 2t a + tj is 1.01 seconds for 
this example. Note that there is also a small contribu- 
tion to the total phase from the adiabatic compression 
and expansion periods, see Fig. [3l The peak density of 
the condensate in its compressed state is 6 x 10 14 /cm 3 
in this case, which is compatible with typical three-body 
loss rates [23|. Shorter gate times could be achieved for 
smaller initial ground state sizes, higher compressed den- 
sities, or larger Feshbach enhancement factors. 

We have shown that a controlled phase of tt be- 
tween individual photons is achievable on the one-second 
timescale under realistic conditions. We hope that our 
proposal will stimulate experimental work in this direc- 
tion. 
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